Optimized Noise Filtration through Dynamical Decoupling 
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One approach to maintaining phase coherence of qubits through dynamical decoupling consists 
of applying a sequence of Hahn spin-echo pulses. Recent studies have shown that, in certain noise 
environments, judicious choice of the delay times between these pulses can greatly improve the 
suppression of phase errors compared to traditional approaches. By enforcing a simple analytical 
condition, we obtain sets of dynamical decoupling sequences that are designed for optimized noise 
filtration and are spectrum-independent up to a single scaling factor set by the coherence time of 
the system. We demonstrate the efficacy of these sequences in suppressing phase errors through 
measurements on a model qubit system, ^Be"'' ions in a Penning trap. Our combined theoretical 
and experimental studies show that in high-frequency-dominated noise environments this approach 
may suppress phase errors orders of magnitude more efliciently than comparable techniques can. 
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Developing improved strategies for maintaining quan- 
tum coherence in the presence of environmental noise is 
central to the advancement of quantum control exper- 
iments and applications. The field of quantum infor- 
mation processing stands to gain significantly from such 
strategies. One approach to maintaining coherence con- 
sists of applying simple pulses intermittently P, 0, S B El 
with computational control operations. Similarly, it is 
possible to design the computationally relevant control 
operations themselves with the same objective in mind 
[a, 0, S]- These techniques are commonly referred to 
as dynamical decoupling. One widely used scheme en- 
tails applying a sequence of Hahn spin-echo-style 7r-pulses 
0, successive pulses being separated by free-precession 
delays. The precise durations of these inter-pulse de- 
lays that would provide optimum error suppression is 
currently a topic of many research efforts [1, 0, SJir 
The Uhrig dynamical decoupling (UDD) sequence 
has shown especially promising improvements over tra- 
ditional approaches such as Carr-Purcell-Meiboom-Gill 
(CPMG) multipulse spin-echo ll| in certain noise envi- 
ronments. 

In this letter we devise and test the noise-suppression 
capabilities of dynamical decoupling pulse sequences tai- 
lored for optimized noise filtration (OFDD sequences). 
To that end, we enforce a simple analytical condition 
that does not depend on detailed knowledge of the noise 
spectrum, yielding, up to a scaling factor in time, a set 
of optimized pulse sequences. Each member of the set 
of sequences is optimized for a unique choice of sequence 
duration and consequently each member has different rel- 
ative pulse locations. In that sense the concomitant dy- 
namical decoupling may be thought of as being "locally 
optimized" This should be contrasted with most 

other spin-echo-style dynamical decoupling, e.g. UDD, 
in which analytical considerations Q naturally lead to 
fixed relative pulse locations for all sequence durations. 

Our studies demonstrate that OFDD sequences sup- 



press errors comparably to, or better than, sequences 
with fixed relative pulse locations in arbitrary noise envi- 
ronments, making orders of magnitude gains in noise en- 
vironments with strong high-frequency components and 
sharp high-frequency cutoffs. The sharpness of the cut- 
off that allows improvement over CPMG by UDD was 
studied in Q. In low-frequency-dominated noise envi- 
ronments OFDD performs similarly to CPMG which is 
an effective slow-noise filter. Finally, we demonstrate the 
efficacy of phase-error suppression by testing OFDD se- 
quences using trapped ions as a model system, showing 
strong agreement between theory and experiment. 



We consider the time evolution of a two-level system 
(qubit) under the influence of random classical noise 
that causes an instantaneous deviation in frequency, P(t), 
from the unperturbed qubit frequency, fl. The dy- 
namics of the qubit are governed by the Hamiltonian: 
H — Ti/2[il. + /3{t)] cTz, where CTz is the Pauli spin opera- 
tor parallel to the quantization axis. Adopt the follow- 
ing conventions for a decoupling sequence of n vr-pulses: 
let T be the total dynamical decoupling sequence du- 
ration, equal to the sum of the free-precession delays 
between pulses, plus the sum of all 7r-pulse durations. 
If the center of the j^^ 7r-pulse occurs at time t = tj 
then let 5j = tj/r and let each 7r-pulse have duration 
Ttt. Assuming the qubit spin- vector is initially aligned 
along the y-axis we use the expectation value of Uy in 
the frame rotating at frequency f2 as a measure of co- 
herence W{t) — \{(Ty{T))\ — e~^^'^\where x(t) is the 
coherence integral, x{r) = f 5(^)^^^du; 
The angle brackets in the former indicate the quantum 
mechanical expectation value and the over-line an ensem- 
ble average. S{ijj) is the power spectral density of (}(t), 
Lo is the angular frequency, and FniyiT) 13, 3] a filter 
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FIG. 1: OFDD sequences for (a) n = 6 and (b) n = 10 pulses. 
The horizontal axis indicates the pulse locations, 5j, relative 
to the total sequence duration r' (vertical axis) that minimize 
Af{t'). The horizontal dashed line at (a) r' — 15.8 and (b) 
t' = 24.2 indicates when the filter function evaluated at the 
cutoff frequency is equal to one. 



function dependent on the n 7r-pulses locations 
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Tailoring the filter function by choosing appropriate Sj 
can lead to enhanced suppression of dephasing. 

The coherence integral is smallest, and hence the er- 
ror suppression strongest, when the overlap between the 
noise power spectrum and the filter function is mini- 
mized. Thus we enforce the condition that the area under 
the filter function Af{t) be minimized over a relevant 
frequency domain [0,W£)] where wd is a potentially un- 
known high-frequency cutoff. This area depends both on 
the relative pulse locations, dj, and the total sequence 
duration r; as such, optimization yields a set of pulse se- 
quences, each sequence optimized for a given value of r. 
In the most general case Ap{t) also depends on Ttt, but 
the following discussion focuses primarily on the instan- 
taneous pulse approximation (t^t = 0), which is appropri- 
ate in qubit systems where 7r-pulse durations axe small 
compared to the system's coherence time, e.g. 
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Transforming to dimensionless units by defining to' = 
w/wd and r' =a;Dr we want yli?(T') = ujd Jq Fn{i^'T')duj' 
to be a minimum. No detailed knowledge of the noise 
spectrum enters into this condition and ud appears only 
as an overall scaling factor, so that the same set of se- 
quences minimizes Ap(t') for any lod. In practice, there- 
fore, the same set of sequences can be used for any noise 
environment by scaling the duration of all sequences in 
the set by the same factor while maintaining the same 
relative pulse locations in each individual sequence. 

For fixed r' we determine the relative pulse locations 
that minimize Ap(t') by running a numerical search rou- 
tine. We use the UDD sequence as an initial guess at 
t' — and the sequence obtained after convergence as 
the initial guess for the next choice of r', and so on. In 



Fig. [T]Ja) and (b) we show the sets of OFDD sequences 
so obtained for a range of sequence durations and assum- 
ing, respectively, n = 6 and n = 10 instantaneous 7r-pulses, 
symmetrically distributed around 0.5r'. 

The observation that the condition of minimized 
Af{t') yields a set of pulse sequences which provide near- 
optimum performance in any noise environment consti- 
tutes the central result of this work. These sequences lead 
to significant gains over existing ones in high-frequency- 
dominated noise environments as we show next. 

We present simulations comparing the error suppres- 
sion capabilities of OFDD sequences obtained by mini- 
mizing Af{t') only, to that of UDD and CPMG, as well 
as the sequences obtained by minimizing the spectrum 
dependent xi''') directly (referred to as LODD - locally 
optimized dynamical decoupling [l2|). In Figs. [21^ a)- (d) 
we plot simulations of the decoherencc, (1 — W{t'))/2, 
having used n = 6 and n= 10 instantaneous vr-pulses, and 
a dimensionless spectrum S{ui')/ujd = culo where a is 
the dimensionless noise strength. We chose two limiting 
cases to illustrate the range of possible outcomes when 
using OFDD sequences. Those are an Ohmic, 7 = 1 (rel- 
evant to a spin-boson model for semiconductor quantum 
dots 0]), and a "1//" (7 = —1) spectrum respectively. 
We assume an infinitely sharp cutoff for the Ohmic, but 
for the 1// spectrum a soft high-frequency cutoff ~ lo'^'^ 
and a sharp low-frequency cutoff at oj' — 0.001 (to pre- 
vent spectrum divergence as to' 0). The main fea- 
tures of this comparison are that, in the high- frequency- 
dominated noise environments. Fig. [H^a) and (b), the 
OFDD approach achieves superior error suppression over 
the entire coherence time. The OFDD error curve has 
roughly the same polynomial growth as that of UDD, 
but extends the coherence time by a factor 1.5. Our 
simulations show that these characteristics persist inde- 
pendently of the noise strength and pulse number up to 
~ 20 TT-pulses, beyond which the optimization becomes 
numerically challenging. Note that the OFDD result dif- 
fers from the LODD approach by a relative error of only 
a few per cent (the two curves are nearly indistinguish- 
able). In the low-frequency-dominated 1// environment. 
Fig- m^c) and (d), all four approaches perform similarly. 
We have verified that similar improvements occur for 
supra-Ohmic noise spectra (7 > 1) with sharp cutoffs, 
results are similar to CPMG for 7 < — 1, while interme- 
diate benefits are achieved for white noise with a sharp 
cutoff. As a rough measure of what constitutes a "sharp" 
cutoff, we find empirically that significant benefits arise 
from our approach if the integrated noise power beyond 
a high-frequency peak is less than ^10% the total in- 
tegrated noise power and contained within a frequency 
width one tenth of lod- 

These numerical results motivate us to develop an 
experimental procedure for implementing OFDD se- 
quences. To that end, consider the following analysis. 
The average relative delay between pulses is sa 1/n. At 
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FIG. 2: Comparison of decoherence, (1 — W{t'))/2, due to 
phase errors when using n = 6 (left) and n — 10 (right) pulse 
CPMG (dotted lines), UDD (dashed lines), OFDD (solid 
lines) and spectrum dependent LODD (dot-dashed lines) in 
an (a), (b) Ohmic and (c), (d) 1// noise environment respec- 
tively. We assumed — and used a = l throughout. 



times t' <C 1 all the phase factors in Eq. ([T]) (within 
the domain Ljj'e[0, 1]) are small compared to tt/2, and the 
different contributions in the sum add destructively (in 
particular i^„(0) = 0). At times r' 1 the phase fac- 
tors in successive terms in the sum can differ by order 
^ TT, and the terms in the sum may add constructively. 
The filter function then oscillates rapidly as a function 
of frequency around its average value An-\-2 and the sys- 
tem is largely dephased. The crossover between the two 
regimes occurs when w/jr/n k, tt, i.e., r' k, nn, leading 
to Fnir') ~ 0{1). This is indicated by the horizontal 
lines in Figs.[TJa) and (b), designating r' = r'p^^, where 
Fnir') = 1. 

With these considerations in mind, appropriate scaling 
of OFDD sequences for a given experimental setting (i.e. 
an effective measurement of lod) can be accomplished ex- 
perimentally using a feedback routine as follows: 

1) Numerically calculate the OFDD set for n-pulses in 
dimensionless units over a time domain [0, 2t^^^]. 

2) Measure the coherence time Tc, i.e. the time for the 
error to increase to 1/e of its asymptotic value, for the 
n-pulse sequence, e.g. using CPMG or UDD. 

3) Associate the pulse sequence (obtained in (1)) at t^^i 
with the measured Tc as a first estimate of the scaling. 

4) Find the sequence in the OFDD set that optimally sup- 
presses errors for t — Tc using a one-dimensional search 
algorithm and experimental feedback. Only a single se- 
quence in the set should provide optimum error suppres- 
sion when the Tc is chosen as the sequence duration since 
the sequences are locally optimized. 



5) Scale the duration of all sequences in the set using the 
result of the above optimization. 

In the above prescription, a single-parameter feedback 
algorithm, implemented at only one choice of sequence 
duration, simultaneously determines the scaling of all se- 
quences within an analytically derived OFDD set. This 
stands in contrast to the n-dimensional feedback opti- 
mization routine (for n pulses) used in our earlier work 
as a means of finding optimized sequences, and which 
must be repeated for each choice of sequence duration. 

The generality of the OFDD approach discussed so far 
applies strictly only to the instantaneous 7r-pulse limit. 
When the pulse durations are nonzero and comparable 
to the inverse cutoff frequency, the relative time scale of 
TTr-to-l/cJc - possibly an experimental unknown - enters 
the problem. However, we find empirically that (up to 
ludTt^ « 1) the pulse sequences obtained by accounting 
for in minimizing Apir') are nearly the same set of 
sequences as shown in Fig. 1, but contracted along the r' 
axis as compared to the instantaneous pulse case. Hence 
the two scenarios differ primarily in the same scaling fac- 
tor that the feedback is designed to determine, and the 
technique still works. After scaling, significant deviation 
between the instantaneous and finite-duration cases oc- 
curs only at short times, when the total free-precession 
time is comparable to uTt^, but there the error suppres- 
sion is extremely strong for any dynamical decoupling 
approach. 

We use ^Be"*" ions in a Penning trap as a model qubit 
system to demonstrate the proposed dynamical decou- 
pling scheme. We present only a brief summary of our 
experimental system, as it has been described in detail 



elsewhere [12|, |13|, [16| . A few hundred to a few thousand 
ions are trapped and Doppler laser-cooled to ~ 1 mK. The 
- 124 GHz (at 4.5 T), 2s^Si/2\mi = 3/2, m,/ = -1/2) = 
i) — > |m/ = 3/2, mj = +1/2) = |t) spin-flip transition of 
^Be"*" serves as a qubit. Coherent rotations of the qubit 
are induced by a microwave field, producing a tunable 
Rabi flopping 7r-pulse duration, here 229 /iS. Each ex- 
periment begins by optically pumping all ions into the 
state (It)) which is bright to cooling light fluorescence. 
The qubits are then rotated to lie along the y-axis by 
applying a 7r/2-pulse before initiating the decoupling se- 
quence. The decoupling sequence ends with a 7r/2-pulse 
that in the absence of dephasing rotates the qubits to 
the dark state (|i)). The accumulation of phase errors 
is manifested as nonzero fluorescence at the end of the 
experiment, due to a nonzero probability of qubit pop- 
ulation in II). The normalized fluorescence count rate 
measured after application of a decoupling sequence is a 
measure of this error and is given by (1 — W{t))/2. We 
achieve a minimum combined operational and measure- 
ment fidelity of ~ 99 % for sequences with rt < 10. A 
desired noise environment is synthesized by frequency- 
modulating the microwave drive as explained in jl3| . 
Figure [3] compares the experimental performance of 
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FIG. 3: Time evolution of the normalized fluorescence counts, 
equivalent to (1 — W{t'))/2, measured using ^Be^ qubits in 
a Penning ion trap in (a) the ambient noise environment and 
a synthesized Ohmic noise environment with (b) a frequency 
cutoff near 500 Hz and (c) a frequency cutoff near 250 Hz. We 
used 6 TT-pulses in each sequence and open circles represent 
CPMG, triangles UDD, and diamonds OFDD. The data are 
well described by the fits (solid lines) which use the analyti- 
cal expression for the coherence, the measured noise spectra 
in both environments, and the noise strength as a free pa- 
rameter. We are at this time unsure of the cause of the minor 
discrepancy between the fits and the UDD and OFDD data in 
(c). This may be the consequence of external noise influences 
not accounted for in the synthesized spectrum. 

different decoupling sequences measured in (a) the ambi- 
ent noise environment of our trap, which approximately 
scales as ~ l/w^ (b) and (c) a synthesized Ohmic 

spectrum with a sharp cutoff around 500 Hz and 250 Hz 
respectively. We used here the OFDD sequences shown in 
Fig.[IJa), having subdivided the interval < r' < 30 into 
about 3000 intervals. After measuring a 6 7r-pulse CPMG 
error curve we choose for each spectrum a fixed sequence 
duration r w Tc {t ~ 20 ms, 2.8 ms and 5.5 ms respec- 
tively) and perform a golden-section search with experi- 
mental feedback to determine which OFDD sequence of 
the set gives optimum error suppression for the selected 
T. The feedback algorithm converges to the optimized 
sequence within ~ 10 iterations. (More specifically to 
a narrow band of sequences all of which give an error 
indistinguishable within the measurement noise. This 
band can in principle be made arbitrarily narrow with 
sufficient averaging during measurement.) This fixes the 
appropriate sequence associated with any other choice 
of T. As anticipated from the simulations presented ear- 
lier, we obtain comparable results for the different decou- 
pling techniques in the ambient noise environment with 
a soft high-frequency cutoff. Fig. EJa). By contrast, in 



the Ohmic noise environment with a sharp cutoff around 
500 Hz, Fig. [3ljb), the OFDD sequence shows significant 
improvements over the other sequences even in the low- 
fidelity regime (decoherence > 1 %). To illustrate that 
the feedback algorithm finds the appropriate scaling for 
arbitrary cutoff we decrease the cutoff by a factor of two 
and repeat the experiment. Fig. [3t^c). This increases the 
coherence time of the system, but the benefit due to 
OFDD remains. Note that despite having finite-duration 
TT-pulses lodTtt ~ 0.7 the technique still works for reasons 
explained earlier. 

In conclusion, we have presented a dynamical decou- 
pling technique designed for optimized noise filtration by 
enforcing an analytical condition that does not rely on 
any detailed knowledge of the noise power spectrum. The 
technique significantly reduces the burden of performing 
measurement-feedback-based optimization of decoupling 
sequences [13], as it relies on only a single parameter 
search at one point of the error curve. Under appropriate 
conditions it improves error suppression by several orders 
of magnitude compared to those of standard sequences. 
Experimental measurements using a model quantum sys- 
tem and artificially engineered noise environments vali- 
date the predicted performance of these sequences. 
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